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Abstract 

We construct solutions to the relativistic Navier-Stokes equations that describe 
the long wavelength collective dynamics of the deconfined plasma phase of = 4 
Yang Mills theory compactified down to d = 3 on a Scherk-Schwarz circle and higher 
dimensional generalisations. Our solutions are stationary, axially symmetric spinning 
balls and rings of plasma. These solutions, which are dual to (yet to be constructed) 
rotating black holes and black rings in Scherk-Schwarz compactified AdSs and AdSe, 
and have properties that are qualitatively similar to those of black holes and black 
rings in flat five dimensional supergravity. 
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1 Introduction 

A particularly interesting entry in the dictionary between gauge theory and gravity links 
deconflned or 'gluon plasma' phase of Yang Mills theory to black branes and black holes 
in gravity. In this paper we study aspects of this connection in the context of speciflc 
examples. In most of this paper we study d = A, SU (N) , A/" = 4 Yang-Mills at 't Hooft 
coupling Qym^ ~ compactifled on a Scherk-Schwarz (the remaining 2 + 1 dimensions 
are non compact). The low energy dynamics of this theory is that of a 2 + 1 dimensional 
Yang-Mills system that undergoes deconflning phase transition at a flnite temperature [1] . 
At large N and strong 't Hooft coupling this system admits supergravity dual description; 
the low temperature conflning phase is dual to a gas of IIB supergravitons on the so called 
AdS soliton background [T] 
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where i = 1, • • • , 2, ^ ~ ^ + 27r, = ancfl 
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T,(m) = 1- -eM . (2) 



The high temperature phase of the same system (at temperature T = is dual to 
the the black brane 



ds^ = L'ol fe^" {-Ti,{u) dt^ + d^2 + d.wf) + du^^ . 



(3) 



The thermodynamics of the high temperature phase are determined in the bulk description 
by the usual constitutive equations of black brane thermodynamics [2] 

For T > Tf. this free energy is negative, and so (in the large N limit) is smaller than the 
(9(1) free energy of the 'confined' gas of gravitons. Consequently, the system undergoes a 
deconfinement phase transition at temperature 7^1^ 

Just as the mean equilibrium properties of the deconfined phase are well described by 
the equations of thermodynamics, the statistically averaged near-equilibrium dynamics of 
this phase is governed by the equations of fluid dynamics - the relativistic generalisation 
of the Navier-Stokes equations. These equations accurately describe the time evolution of 
fluid configurations whose space time derivatives are all small in units of the mean free 
path, which is of the same order as the mass gap of the theory [2l[3]. The same equations, 
augmented by appropriate surface terms, may also be used to study the dynamics of large 
lumps of plasma localised in the gauge theory vacuum. 

The properties of the surface that separates the plasma from the vacuum, may be 
studied in the context of the simplest plasma profile with a surface; a configuration in 
which half of space, x < 0, is filled with the plasma. The surface at x = is a domain wall 
that separates the plasma from the vacuum. The net force on this domain wall vanishes 
(and so the system is in equilibrium) when the plasma that fills x < has vanishing 
pressure, i.e. at T = 7^ in the large limit. The bulk gravity dual of this solution was 
constructed numerically in [2]; this configuration interpolates between the black brane at 
T = % for X < and the vacuum at x > 0, via a domain wall. The thickness and surface 
tension of this domain wall may be read off from this gravitational solutions, and were 
estimated, in [2] at approximately 6 x -^^^ and a = 2.0 x '^'^2'^" ■ 

More generally, one would expect a finite lump of plasma that evolves according to 
the relativistic Navier-Stokes equations map in the bulk to a 'black hole' that evolves 
according to the Einstein equations. Provided all length scales in the plasma solution are 



^Notice that, at large u, Tx{u) ~ 1, so ([T]) reduces to AdSd+2 in Poincarc-patch coordinates, with u as 
the radial scale coordinate, and with one of the spatial boundary coordinates, 0, compactified on a circle 
(the remaining boundary coordinates, Wi and remain non-compact). 

'^Tc = l/27r in the dimensionless imits of ([3]) 
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small compared to the gauge theory mass gap (which is of the same order as the domain 
wall thickness), the dual bulk solution is well approximated by a superposition of patches 
of the black brane solution (with temperature varying across the patches) in the bulk and 
patches of the domain wall solution described in the previous paragraph. It follows (at least 
for stationary solutions) that the 3 dimensional black hole horizon topology (at any given 
time) is given by an (physically this is the 6 circle) fibred over the two dimensional fluid 
configuration at the same time, subject to the condition that the contracts at all fluid 
boundaries. Consequently, fluid configurations with different topologies yield bulk dual 
black hole configurations with distinct horizon topologies. We will return to this point 
below. 

This paper is devoted to a detailed study of certain 'stationary' configurations of the 
plasma fluid; i.e. time independent, steady state solutions to the relativistic Navier-Stokes 
equations. The simplest configurations of this sort was studied already in [2]; the plas- 
maball is a static, spherically symmetric lump of fluid at constant local pressure P with 
P = a/R where R is the radius of the lump and a its surface tension. In this paper we 
study the more intricate spinning lumps of stationary fluid. These lumps carry angular 
momentum in addition to their mass. 

It turns out that the relativistic Navier-Stokes equations admit two distinct classes of 
solutions of these sort. The first class of solution is a simple deformation of the static plas- 
maball; it is given by plasmaballs that spin at a constant angular velocity. The centripetal 
force needed to keep the configuration rotating in this solution is provided by a pressure 
gradient. The local plasma pressure (and hence local temperature and density) decreases 
from the edge (where it is a positive number set by the radius, surface tension and rotation 
speed) to the centre. As large enough angular velocity the pressure goes sufficiently nega- 
tive in the core of the solution to allow for a second kind of solution of these equations; an 
annulus of plasma fluid rotating at constant angular velocity u. The local plasma pressure 
is positive on the outer surface and negative at the inner surface; the numerical value of 
the pressure in each case precisely balances the surface tensions at these boundaries. 

We now describe the moduli space of spinning plasmaball and plasma ring solutions 
in a little more detail. In figll](a) we have plotted the energy-angular momentum plane, 
which we have divided up into 4 regions. In region C (low angular momentum at fixed 
energy) the only rigidly rotating solution to the equations of fiuid dynamics is the rotating 
plasmaball. At higher angular momentum (region B) in addition to the rotating plasmaball 
there exist two new annulus type solutions which we call large and small ring solutions. As 
their names makes clear, the solutions are distinguished by their size; the large ring has a 
larger outer radius than the small one. On further raising angular momentum (region A) , 
the small ring and the ball cease to exist; in this region the large ring is the only solution. 
Finally, at still larger angular momentum (region O) there exist no solutions. 

In fig{T](b) we have plotted the entropy of the three different kinds of solutions as a 
function of their angular momentum at a particular fixed energy. At angular momenta for 
which all three solutions coexist (region B) the entropy of the small ring is always smaller 
than the entropy of either the large ring or the black hole. Upon raising the angular 
momentum, the solution with dominant entropy switches from being the ball to the large 
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O - no solutions, A - large ring, B - large ring, small ring and ball, C - ball. 




Figure 1: (a) Regions where ball and ring solutions exist, (b) their entropy as a function 
of angular momentum at fixed energy. 



ring; the first order transition between these solutions occurs at an angular momentum 
that lies on a 'phase transition line' in the bulk of region B. This picture suggests - and we 
conjecture - that the ball and the large ring are locally stable with respect to axisymmetric 
fluctuations, while the small ring is locally unstable to such fluctuationslfl In § §5.51 we 
perform a 'turning point' analysis of our solutions, to find some evidence for this guess. 

Let us now turn to the bulk dual interpretation of our solutions. The fluid for the 
spinning plasmaball is topologically a disk; consequently the horizon topology for the dual 
bulk solution - the 5"^ fibration over this disk - yields an 5''^. The bulk dual of the spinning 
plasmaball is simply a rotating five dimensional black hole. On the other hand the fluid 
configuration of the plasmaring has the topology of x interval; the fibration over this 
configuration yields x S"^; i.e. a five dimensional black ring! Notice that in addition to 
the isometry along the S^, these ring solutions all have a isometry on the S"^ corresponding 
to translations along the Scherk-Schwarz circle. This additional isometry, that does not 
appear to be required by symmetry considerations, appears to be a feature of all known 
black ring solutions in flat space as well. 

Using the gauge theory / gravity duality, the quantitative versions of the figH] give 
precise quantitative predictions for the existence, thermodynamic properties and stability 
of sufficiently big black holes and black rings in Scherk-Schwarz compactified AdSs spaces. 
While these gravitational solutions have not yet been constructed, their analogues in flat 
5 dimensional space are known, and have been well studied. The general qualitative fea- 
tures (and some quantitative features) of figHJare in remarkably good agreement with the 
analogous plots for black holes and black rings in flat five dimensional space (see § §5.41 for 
a detailed discussion). 

The constructions we have described above admit simple generalisations to plasma 
solutions dual to black holes and black rings in Scherk-Schwarz compactified AdSg spaced 

■^It is possible that the large ring exhibits Plateau-Rayleigh type instabilities that break rotational 
invariance; such modes would map to Grcgory-Laflammc type instabilities of the bulk solution (see also [4]). 
We thank T. Wiseman for suggesting this possibility. 

^Note that the spinning plasmaring has no analogue in 1+1 dimensional fluid dynamics, for the excellent 
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Figure 2: Spinning ball and ring solutions. 



As the qualitative nature of the moduli space of black hole like solutions in six dimensional 
gravity is poorly understood, this study is of interest. The boundary duals of these objects, 
in the long wavelength limit, are stationary solutions to the equations of fluid dynamics 
of a 4 dimensional field theory. In ^ we construct such solutions. It turns out that these 
solutions occur in two qualitatively distinct classes. The simplest solutions are simply 
spinning balls of plasma; the fact that these balls spin causes them to flatten out near the 
'poles'. As these balls are spun up, their profile begins to 'dip' near the poles (see figE]). 
As these balls are further spun up, they pinch off at the centre and turn into doughnut 
shaped rings (see figl2D- 

As in the three dimensional case, the horizon topology of the black objects dual to the 
rotating plasmaballs and plasmarings described above, is obtained by fibering the fluid 
configuration with an that shrinks to zero at the fluid edges. This procedure yields a 
horizon topology S'^ for the dual to the rotating plasmaball, and topology x for the 
dual to the plasmaring. As plasmaball and plasmaring configurations appear to exhaust 
the set of stationary fluid solutions to the equations of fluid dynamics, it follows that 
arbitrarily large stationary black objects in Scherk-Schwarz compactified AdSe all have 
one of these two horizon topologies. x is an example of another topology one could 
have imagined for black objects in this space; these would have been dual to hollow shells 
of rotating fluid; however, there are no such stationary solutions to the equations of fluid 
dynamics. 

The analysis of four dimensional fluid configurations, described above, demonstrates 
the power of the fluid dynamical method. In simple contexts, the Navier-Stokes equations 
are much easier to solve than the full set of Einstein's equations, and rather easily reveal 
interesting and nontrivial information. It would be interesting to extend our analysis of 
fluid dynamical models in various directions to obtain information about the moduli space 
and stability of classes of black solutions in AdS spaces. An obvious extension would be 
to move to higher dimensions. As a first step in this direction, we have obtained and 
partially solved the fluid flow equations in 5 dimensional spaces. A complete analysis of 

reason that there is no spin. This tahics with the fact that there are no black rings in four dimensions (at 
least in flat space). 
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these equations would yield the spectrum of black holes in Scherk-Schwarz compactified 
AdSy spaces, in terms of the fluid dynamics of the deconfined phase of the M5 brane theory 
on a Scherk-Schwarz circle. 

Finally, we should point out that there has been a long history within the General 
Relativity literature of treating black hole horizons as surfaces associated with fluids. In 
one of the most recent discussion within this framework, the authors of [5] have modelled 
spinning black holes in d + 1 dimensional (fiat space) gravity by + 1 dimensional lumps of 
incompressible fluid; here the fluid surface represents the black hole horizon. Within this 
framework the 4+1 dimensional black ring, for instance, is modelled by a 4+1 dimensional 
stationary fluid lump of topology x [6]. This description is rather different from 
the AdS/CFT induced description of black rings in Scherk Schwarz compactified AdSs 
as a 2+1 dimensional annulus of fluid. It would be interesting to better understand the 
interconnections between these approaches. 



2 Fluid mechanics and thermodynamics 

In this paper we study aspects of the dynamics of the deconfined plasma described in 
the previous section. A full accounting for the dynamics of the 'gluon plasma' is a very 
complicated problem. However, when the thermodynamic potentials and velocities vary 
over length scales large compared to the quasiparticle ('gluon') mean free path admit 
an effective description in terms of the equations of fluid dynamics. The variables in this 
description are the local values of the plasma or fluid velocity u^{x) and the plasma density 
p{x) (the equation of state, as discussed in § §2. 41 may be used to trade the density for the 
pressure or the temperature). The equations of fluid dynamics are simply a statement of 
the conservation of the stress tensor 

S/^T^"" = df.T"" + T'l^xT^'' + T^^xT^^ = • (5) 

All input of the dynamical nature of the fluid that undergoes this flow appears in the 
specification of the stress tensor in terms of the velocity and density of the fluid, and from 
the thermodynamic equation of state (which determines the pressure and temperature as 
a function of density). In the rest of this brief subsection we describe the functional form 
of the fluid dynamical stress tensor for our plasma fluid in detail. 
The stress tensor can be split into three parts: 

perfect ' dissipative ' surface" 

The first part, T^^^^^^t^ stress tensor for a perfect fluid with no dissipative forces. It 

is a function only of fluid velocity and thermodynamic quantities in the rest frame, and 
not of their space time derivatives. 

The second part, T'dissipativc receives contributions from viscosity and heat flow. In 
the long wavelength limit this piece is linear in the first derivatives of the velocity and 
temperature. 
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The third part, T^^^^.^^^, represents surface contributions to the stress tensor, and re- 
quires more explanation. Any fluid configuration with a surface has large variations in (for 
instance) the fluid density over the scale of the mean free path, in directions normal to the 
surface (see for instance [2]). As a consequence it is impermissable to use the Navier-Stokes 
equations for the fluid in the neighbourhood of the surface. When the deviations of the 
surface from a straight line are small over length scale of the mean free path, however, all 
effects of the surface may approximately be captured by a delta function localised 'surface 
tension' contribution, T^^^^^^c^ to the stress tensor. In the long wavelength approximation, 
this term depends only on the gradients of the surface and not its curvature. 



2.1 Perfect fluid stress tensor 

The most general ultralocal stress tensor one can build out of the fluid velocity and ther- 
modynamic quantities that reduces to what is expected for a fluid at rest is [3, ch.22] 

^pTrfect = (P + P)^'^" + ■ (6) 



2.2 Dissipative part 

Realistic fluids have a dissipative component to their stress tensor in addition to the perfect 
fluid piece. In the long wavelength limit, this stress tensor is a function of the accelera- 
tion, expansion, projection, and shear tensors (see e.g. [3 Exercise 22.6-7] and references 
therein). 



= vx, 

pf.u ^ g^,u ^ (7) 

2 ^ ' a— 1 

In terms of these quantities and the heat flow vector (see immediately below) 

^d^iLipativc = -ceP'"" - '^m'" + (1'^" + ^'<f , (8) 

where C, is the bulk viscosity, rj is the shear viscosity. The heat flux vector, 

g/^ = -KP'^-(5,T + a,T). (9) 

is the relativistic generalisation of g = — kVT (here k is the thermal conductivity); the 
extra term in (jHl) is related to the inertia of flowing heat. 



2.3 Surface contribution 

We will use a simple model of surface tension where the energy stored in the surface and 
the force per unit length are both given by cr, which we take to be the surface tension 
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at the critical temperature computed in [2]. We will ignore any dependance o could have 
on the fluid temperature. This approximation is valid when the fluid temperature at the 
surface does not deviate substantially from 7^. 

Consider a localised lump of fluid whose surface in space is given by the equation 
]{x) = 0. The surface contribution to the stress tensor will be proportional to (T6{f). In 
the long wavelength limit it will only depend on the flrst derivatives of /. The most general 
stress tensor we can build is 

^surface = [« ^7 + /? U^^U^ + 7 {u^dj + d^^fu^ + 6 9^ ^Kf) ■ 

As = — 1 and u^^d^f = (the surface moves with the fluid), the only invariant quantity 
that 7,5 can depend on is (9^/(9^/). We can flx this dependence by demanding 
invariance under reparameterisations of the surfac^ (e.g. f{x) — g{x)f{x), so that df 
gdf + fdg = gdf at the surface). Deflning = J^j^ '- 



^surface = + Bu^u^ + C {u^ T + Fu^ + Dg^^] a^/dfdf5{f) . 

We can flx A, B, C, D by looking at a fluid at rest, u'^ = (1, 0, 0, . . .), with a surface f{x) = x 

fB-D C 0' 

rs"urface= C A + D 0\a5{x)=\0 \a5{x). 
\ 

This gives 



2.4 Equations of state 

To solve the equations of fluid mechanics, one also needs expressions for the various coeffi- 
cients that appear in the stress tensor above in terms of the density. For our purposes, we 
only need to know the thermodynamic properties of the fluid, which could be determined 
from the static black brane solution In this subsection we discuss the free energy, 
temperature etc. of the plasma at rest. This is different from the free energy, temperature 
etc. of the plasmaball/plasmaring. 

For a conformal theory in d dimensions with no conserved charges, dimensional analysis 
and extensivity determine 

= -aVr\ (11) 

with a an arbitrary constant. In our situation, the plasma is dual to the same black brane, 
so it doesn't know about any capping off in the IR except that the energy is measured with 
respect to a different zero. Before reducing on the Scherk-Schwarz circle, it behaves like 

■'However, we should choose a parametrisation such that d^f is weh behaved at the surface, e.g. f = x, 
but not f ~ or f — \fx. 
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a conformal theory in d + 1 dimensions plus a vacuum energy density. After dimensional 
reductioij^, we have 

T = V[p^- al"^^) . (12) 

This gives 



S = J^ + TS = V{po + daT'^^^) . 
In terms of intensive quantities, we have 

P = ^ , P + P=[-^){p-Po). 

or, in three dimensions 

P=^, P + P = ^(P-P0), 

4aV^ .3/4 ^_^p-po^^/^ ^''^ 



^- 33/4 (P-Po) , V 3« 

Note that the critical density and temperature are those for which the pressure is zero 

i/(d+i) 



p,= (d+l)po, '^c=(^) . (16) 



For the black-brane equation of state 



However, the values of these constants will not be important below. 



3 Rigidly rotating configurations 

In this section, we study stationary, axially symmetric rotating fluid configurations, whose 
equation of state is presented in various forms in § §2.4[ We choose the axis of rotation as our 

^Strictly speaking, it is not a dimensional reduction as we will have plasma temperature of the same 
order as the Kaluza-Klein scale. Rather, we are restricting attention to classical solutions that do not vary 
in this compact dimension. 
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origin in polar coordinates; in these coordinates the fluid density is a function only of the 
radial coordinate r, and the (t, r, 0) components of the velocity are given by = 7(1, 0, uj) 
with 7 = (1 — cj^r^) We will find two distinct kinds of solutions; rotating plasmaballs 
with the topology of a two dimensional disk, and plasmarings with the topology of a two 
dimensional annulus. The configurations we find are exact solutions to the equations of 
relativistic fluid dynamics; in § §5.21 we will demonstrate that these equations accurately 
represent plasma dynamics for large enough plasmaballs and plasmarings. 

3.1 Equations of motion 

Our fluid propagates in flat 2+1 dimensional space. In polar coordinates 

ds^ = -dt^ + dr^ + r^d^^ . (18) 
This gives the following non-zero Christoffel symbols: 

r;, = -r rt, = rj, = i . (19) 

For the stationary, axially symmetric configurations under consideration, dtT^^ = 
d^Ti''' = 0. Using (dnD, ® becomes 

= VnT"' = drT' + -T\ (20) 

r 

= V.T'^^ = 9^T^^ + -T'"' - rT''"^, (21) 

r 

= V.T^'^ = a^T""^ + -T""^. (22) 

r 

The boundaries are /„ = r — r„, with n labelling the different boundaries (outer for the 
disk, outer and inner for the annulus). The 'perfect fluid part' of the stress tensor is 

h\p + uj^r^P) ^Mp + P) \ 
V iMp + P) l^(^Vp + P)/ 

and the surface stress tensor 

/I \ 
\0 -i/ 

For the dissipative part of the stress tensor, we find 6 = a^" = and 



d,r + a,r= ( 0,7^ 
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,0 , (25) 
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so that 

dissipativc ' 



'T' 




1 






■ 





°) 


.7. 


Vo 


UJ 


0/ 



(26) 



We will now write the equations of motion V^T^*^ = temporarily ignoring the contri- 
bution from this heat flow, T'^ssipativc! ^^^^ turn out (we see this immediately below) that 
^dissipative actually vanishes on our solutions, justifying this procedure. 

The only non-trivial equation of motion, (12T|) . can be written as 

^ = r^(p + ^)-E-^(---^)- (27) 

dr 1 — cu^r^ r 

For these fluids (with no chemical potentials for any conserved charges), P = —f is a 
function only of T, so P + p = sT and ^ = s- So, away from the boundaries, ( 127|) 
becomes 



d 
dr 



^ (28) 







It follows that T'jjggjpj^^.j^g vanishes for rigid rotation, justifying our neglect of heat flow. 

Our discussion has not assumed a specific form of the equation of state. Using this 
particular equation of state of our plasma ( ITSl) . we can rewrite (l28l) in the fluid interior as 

(p(r) - po) (1 - uj^r^Y = constant. (29) 
Integrating (1271) across a surface gives 

P> - P< = -- . (30) 
r 

where P> and P< are the pressures at infinitesimally greater and smaller radii than the 
location of the surface. 

3.2 Spinning ball 

Let us first study a fluid configuration with a single outer surface at r = Tq with P> = 0. 
Using the equation of state ffT^ . the boundary condition flHUl) can be written as 

p(ro) = 4po + — . (31) 

If we define dimensionless variables 

UJ = — , r = , V = ujr = ur , (32) 

Po cr 
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then (1291) can be written as 



Note that the range of v is [0, 1] and p{v) — po is always positive for this solution, as is 
required for the last equation of ( fT5i) to make sense. 

We can also compute the local plasma temperature using ( |T5l) 



r = 7(^^^)'''. (34) 

3.3 Spinning ring 

We now turn to solutions that have an inner surface well an outer surface. In addition to 
the boundary condition at the outer radius flHTl) we now have 

p{n) = 4po - - . (35) 



So the full solution is 



3po 



(36) 



Note that p{v) — po > provided that v\ > lj. 

Again, we can compute the local plasma temperature using (fT5|) 

r = 7 (^M±M^ = ^ i^MzM^ . (37) 

The two functions, g±{v) are schematically plotted in figOlfor some value of u, where 
we have labelled special velocities v* and v*. As Vi < 1, it is necessary that u; < 1. 

We can see that there are no solutions to g+{vo) = g-{vi) for Vo < v* and two solutions 
for Vo > V*. One of these has Vi < v* (the thick ring) and one has Vi > v* (the thin ring). 
The distinction between 'thin' and 'thick' rings will not prove physically important. In 
§ §5.11 we will find it physically useful to distinguish between distinct ring solutions (we 
will call these large and small rings) at the same values of conserved charges (energy and 
angular momentum), rather than the parameters Vq and uj. 
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9-{v) 



1 



Figure 3: Graph of g±{v) showing possible values of v, 



4 Thermodynamic Potentials 

In this section, we compute the thermodynamic potentials (energy, angular momentum, 
entropy, etc.) for the spinning plasmaball and plasmarings themselves, rather than their 
constituent plasma. This includes contributions from the kinetic energy of the plasma as 
well as its internal energy. 

The constitutive relations we find are predictions for, e.g., entropy as a function of mass 
and angular momentum of the dual gravity solutions. 

4.1 Densities 

In this subsection we list formulae for energy density, angular momentum density and 
entropy density. In the next subsection we will integrate these expressions to find explicit 
formulae for the energy, angular momentum and entropy of spinning plasmaballs and 
plasmarings. 

The energy density is given by 



T" = 7' (P + ^V^P) + ^ a 5(r - r„) 



n 





+ ^225t;5(w2 




n 




(39) 



The entropy density is given by 




(40) 
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4.2 Integrals 

We can define some dimensionless variables 



S 



5/4 e 



T = T\ — 



1/4 



f2=^. (41) 
Po 



The last of tfiese ensure tliat T 

4.2.1 Spinning ball 

Energy 



i)z 1^ = (f follow from T = (f )^ and 



— + 5a;fo — 



E 



Angular momentum 
Entropy 



4.2.2 Spinning ring 

Energy 



E 



Angular momentum 



A{vl - ) - {vl - vf) + huj{vo + v{) - uj{vl + ) 



Entropy 



5^4 



l-vl\l + - 

Vo 



uj--' 



3/4 



1 - 1 - - 



3/4- 



4.3 Temperature and angular velocity 



(42) 
(43) 
(44) 



(45) 



(46) 



(47) 



In this subsection we determine the temperature and angular velocity of spinning plasma- 
balls and plasmarings using 



T - I — 



dE 



(48) 
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Figure 4: Scatter plots of energy and angular momentum of (a) spinning balls and (b) 
rings. 



Note that the temperature defined above is different from the local plasma temperature, 
T (which varies across our solutions), in fll 311151) : the angular velocity defined above will 
turn out to be uj on all our solutions, although it is apparently a priori different. 

It may be verified that the temperature and angular velocity of plasmaballs is given by 

f=[g+{v,)Y'\ n = u. (49) 
The corresponding expressions for the rings are identical 

f = [g+MY/' = [gM)?^\ ^ = ^- (50) 

Thus, local temperatures, %, and angular velocities, u, for both the ball and the ring, 
are given simply in terms of T and Q 

T 

T = , , UJ = Q . 

5 Solutions at fixed energy and angular momentum 
5.1 Existence 

In figJH we display scatter plots for the energy and angular momentum of ball and ring 
solutions over the full range of solution parameters|l| The various regions of existence 
of the plasmaball, thin plasmaring and thick plasmaring in the E — L plane are drawn 
schematically in figJS] 

The ball solution exists over a region C in the E — L plane. At the boundary of the 
region C the ball solution attains its maximum value of unity. Using (14211431) we find an 
analytic expression for the boundary of C: 

L = — {3E + 4)3/2 _ 9^ _ 8 _ for large E . (51) 

27 L J 3'^'^ 



^In order to generate these plots for the baU, a range of vales of {vo,i^) were chosen and {E,L) were 
computed using (|42l43p . For the ring, a range of vales of (wi,w) were chosen, Vq was computed using ([55)1 . 
and [E,L) were computed using (|45l46p . 
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O - no solutions, A - 1 ring, B - 2 rings, C - 1 ball. 




E 



Figure 5: Regions where (a) ball and (b) ring solutions exist. 



From poj) . we see that balls on this boundary saturate the extremality bound (i.e. have 
zero temperature). 

Like the balls, rings of a fixed energy have a maximum value of angular momentum. 
Rings at the edge of this bound (the boundary between O and A in fig|5b) have Vq = v\ = 1 
and so are extremal (see ( l50l) ) and of zero width. Using (145 f 46 1) the 0-A boundary is given 
by 

i = ^. (52) 

(this expression is valid only for > 8, L > 2; at lower energies u exceeds unity). 

As we lower angular momentum of the solution, this ring moves away from extremality 
and increases in width. At a particular angular momentum (the boundary between region 
A and region B) a new ring solution comes into existence. The corresponding solution has 
t'o = 1, v\ = and so is extremal (see ( l50l) . Using fl45|46p . the analytic expression for the 
A-B boundary is given by 



^=27 



(3E + if'^ - + ij ~ for large E , (53) 



(for > 8, L > 2 as above). In the high energy limit E ^ 1 the ratio of angular 
momentum for the new extremal rings (at the A-B boundary) and extremal plasmaball 
tends to unity, (even though the the difference between angular momenta does not go to 
zero). Consequently the leading high energy behaviour of (1531) and (l^T]) is the same in this 
limit, as is also clear from figH We emphasise that, for our solutions, the extremal ball 
and extremal thick ring are not quite identical (as is the case for black holes and small 
black rings [8] in fiat space) as the inner radius of our extremal thick rings does not vanish. 
However, the inner radius of the extremal thick ring is always (for all values of energy) of 
the same order as the thickness of the domain wall. As the fiuid dynamics approximations 
fail precisely under these conditions, it could well be that the new extremal plasmaring 
and extremal plasmaball are actually identical configurations. 
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70 SO ^90 

L 




Figure 7: Regions where ball and ring solutions exist. 



As we further lower angular momentum, the new ring solution moves away from ex- 
tremality; this new solution always has a smaller outer radius than the 'original' ring 
solution (the solutions that also exists in region A), as shown in figEl As a consequence, 
we refer to these two ring solutions as small and large respectively. 

Further lowering angular momentum, we hit the boundary between regions B and O 
where the two ring solutions merge into each other. At still lower angular momentum, we 
have no ring solutions. The C-0 boundary may thus be obtained by minimising L at fixed 
E. If one uses ( 136|) to eliminate u5, this amounts to 




The existence of plasmaball and plasmaring solutions in the E-L plane may thus be 
summarised as in figE 

5.2 Validity 

As we have described above, plasmaballs and plasmarings are exact solutions to the rel- 
ativistic Navier-Stokes equations (supplemented by sharp surface boundary conditions). 
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However these equations of fluid dynamics accurately capture the dynamics of the fluid 
plasma only under certain conditions. In our discussion we have assigned a well defined 
pressure and temperature to the fluid at each point in space. Clearly this procedure is 
valid only when the variation of these thermodynamic quantities is small over the length 
scale of the mean free path of the quasiparticles (roughly gluons) of our system. The mean 
free path is of the same order as the mass gap of the theory, which in turn is similar to the 
deconfinement temperature 

Second, we have treated the surface of the plasma as sharp; in reality this surface has a 
thickness of order 7^~^. Consequently, our treatment of the surface is valid only when its 
deviation from a straight line occurs on scales large compared to 7^~^ (higher derivative 
contributions to the surface stress tensor, which we have ignored in our treatment, would 
become important if this were not the case); further we must also require that only a small 
fraction of the fluid should reside in surfaces. 

Thirdly, we have ignored the fact that the surface tension is a function of the fluid 
temperature at the surface, and simply set a = ct(7^). This is valid provided that T /% ~ 1 
at all surfaces. 

Finally, the fluid evolution equations, by their very nature, track mean velocities and 
energy densities, ignoring fluctuations. In our context this approximation is justified by 
large A^; fluctuations are suppressed by powers of l/N"^, dual to the suppression of quantum 
metric fluctuations in the bulk. 

Recall that 

Pq^N^%\ ar^N^T,^ so — ~ — 

% Po 

(for the domain wall solution of [2], — = 2.0 x ^ and the thickness is 6 x tt^). 

As an estimate of the scale over which thermodynamic quantities vary, we compute the 
fractional change in the fluid temperature over the distance of a mean free path. As the 
temperature is proportional to 7, we should look at 

Id a iJ^r OJV 

— — ln7 ~ 



Tc&r Po 1 — cu^r^ 1 — 

As this takes its maximum value at the outer surface, the condition for the validity of the 
equations of fluid dynamics may be estimated to be 

A«^-^«l. (54) 

Our treatment of the surface as a zero-thickness object is valid if 

1 a 



{ro,ri,ro - r;} > 

Po 



^It was argued in [2] that the mean free path does not scale with N. According to quasiparticle kinetic 
theory, the mean free path is approximately the ratio of the shear viscosity and the energy density. The 
computations reviewed in [3] show that this quantity is of order 1 in the limit of the 't Hooft coupling 
A —>■ 00. 
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(for the ring, the Tq inequahty in the equation above follows automatically from the either 
of the other two inequalities). This condition can be rewritten in terms of our dimensionless 
variables as 



Wo > Co" or 
{fi, Vo — Vi} > or 



To ^ 1 for the ball, 
{^i; — f\} ^ 1 for the ring. 



(55) 



In figJHl we have plotted ln(l/AM), In(fi), ln(fo — ri) and In(fo) for the thin ring, thick 
ring and ball. From the figure we observe that these quantities are large (and so the fluid 
dynamics approximations of this paper are accurate) when our rings and balls have large 
energy and we stay away from the extremality bounds. 

Finally, validity of our approximation of the surface tension as a constant (independent 
of temperature) requires that the maximum and minimum values of ln(T/7^) (which occur 
at the outer and inner surfaces respectively) are both small. We have plotted these quan- 
tities in figM It is clear from these figures that this condition is fulfilled for large energy 
and angular momentum provided that we are not near extremality. 

5.3 Global stability and phase diagram 

Recall that (see figJTj) at fixed values of energy and angular momentum, we have either 0, 
1 or 3 plasmaball / plasmaring solutions. At those values of charges for which multiple 
solutions exist, it is natural to inquire which of these solutions is entropically favoured. In 
fjgJTOVa) we have plotted the entropy of plasma ball and plasmaring solutions as a function 
of angular momentum at fixed energy. 

Note that, when it exists, the small ring always carries lower entropy than both the 
big ring and the plasmaball. At low enough angular momentum the plasmaball is the only 
solution. This solution continues to be entropically dominant (upon raising the angular 
momentum) over an interval, even after the new ring solutions are nucleated. At a critical 
angular momentum, however, the entropy of the large ring equals and then exceeds the 
entropy of the plasmaball (all three ring solutions continue to exist in a neighbourhood 
about this point). The large ring is the entropically dominant solution at all larger angular 
momenta. 

The phase boundary can be seen in figfTUlfb). 

5.4 Comparison with black rings in flat 5D space 

As we have explained in the introduction, the plasmaball and plasmaring solutions of this 
paper are dual to black holes and black rings in the background ([T]). Unfortunately the 
corresponding gravitational solutions have not yet been constructed; however exact Black 
ring solutions to the vacuum Einstein equations in 5 dimensions, were obtained in [8] (see [9] 
for a review). These solutions were further studied in [10]. In this subsection we compare 
the properties these black rings and black holes with our plasmaballs and plasmarings, and 
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Figure 8: Plots of i) ln(l/An) ii) In(ri), iii) ln(ro — r[), iv) In(ro) for (a) large rings, (b) 
small rings and (c) balls. 




Figure 9: Plots of (a) ln(7^j„/7^) (b) \n(Tmax/%) for i) balls, ii) small rings, iii) large rings. 
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Green - ball, Red - large ring, Blue - small ring. 




Figure 10: (a) Entropy, S, as a function of angular momentum, L for fixed energy, E = 40. 
(b) Phase boundary with existence boundaries. 

find broad qualitative agreement between the twoj^ 

In figlH] we have presented a schematic plot for the existence of black hole and black 
ring solutions in 5 dimensional fiat space. This figure looks fairly similar to figs JTllTOl fb). 
The major qualitative difference is the absence of the analogue of the region O (see fig. [7]) 
in fig. [TTJ Thus unlike thin black rings in fiat 5 dimensional space, plasmarings (and so 
black rings in Scherk-Schwarz AdSs) have an upper bound to their angular momentum at 
fixed energy0 

It is interesting to pursue the comparison between these solutions in more detail. The 
gravitational analogue of figHO] (presented as [HI fig-3]) looks fairly similar to our figure. 
The main qualitative differences are: unlike for plasmarings, the entropy of the large fiat 
space black ring doesn't go to zero at a finite angular momentum (it asymptotes to zero 
at infinity) and the entropy of the small fiat space black ring and black hole go to zero at 
exactly the same point instead of the slightly different values that we see. We expect that 
first of these differences reflects a physical difference between black rings in flat space and 
Scherk-Schwarz compactified AdSs, the second difference is an artefact of the breakdown 
of the fluid dynamics approximation for extremal small rings (whose inner radius is always 
of order the mean free path). 

In even greater detail, we could quantitatively compare the boundaries between re- 
gions O, A, B and C (see figIS]). These curves, as well as the phase boundary, may be 
parameterised by L = xE^ at large energies. 

For black holes and black rings in flat space i/ab = Hoc = Vbo = U phase = f • For our 
plasmaballs and plasmarings, as one can see in figiI2] (or from 05111531) for the first three), 
for large energy, we get yoA = 2, i/ab = Hoc = |, Vbo = 1-25 and i/phase = 1-25 (see table 

^While we expect the properties of plasmaballs and plasmarings to match quantitatively with those 
of black holes and black rings in the background ([T|), we could not hope to find better than qualitative 
agreement with the properties of the same objects in flat space. 

^°This upper bound was expected for black rings in AdS. The negative cosmological constant has a 
similar effect to the dipole charge of [TT]. We thank R. Emparan for explaining this to us. 
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A' - thin black ring, B' - thin black ring, thick black ring and black hole, C - black hole. 




Figure 11: Existence regions and phase boundary for black holes / rings. 




Figure 12: (a) log-log plots of the boundaries, (b) gradients of log-log plots, (c) L/ —>■ x. 

It is meaningless to compare the x's directly, as they are dimensionful quantities. How- 
ever, when two |/'s have the same value, the ratio of the corresponding x's is dimen- 
sionless and may be compared. For black rings xab = v/32G/277r, xbo = a/G/tt and 
Xphase = ^256G/2437r, so xab/xqc = 1, xqc/xbo = v^32/27 and xso/xphase = 9^3/16. 
For plasmaballs and plasmarings, if we used the dimensionless quantities ( HTl) . we find 
xoA = jQ, Xab = xoc = 2/3^/^ xbo ~ 0.60 and Xphase ~ 0.67. Therefore xab/xqc = 1 
and Xbo/ Xphase -0.91. 

This is summarised in table [TJ Note that the extremality boundaries, OA, AB and OC, 
occur precisely where at least one of the approximations discussed in § §5.21 breaks down. 
Therefore, nothing quantitative about these boundaries should be trusted. 

5.5 Turning point stability 

We have seen in § §5.31 that the spinning plasma solution of maximal entropy is the plas- 
maball (at low angular momentum) or the large plasmaring (at high angular momentum). 
The 'phase transition' between these two solutions may be thought of as being of first order 
(in the sense that the two competing solutions are different at the phase transition point). 
The small plasmaring is entropically subdominant to both the plasmaball and the large 
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Quantity 


Black rings 


Plasmarings 


yoA 


N/A 


2 


Vab 


3/2 


3/2 


yoc 


3/2 


3/2 


Vbo 


3/2 


1.25 


Vphase 


3/2 


1.25 


xab/xoc 


1 


1 




v/32/27 


N/A 


■^BO / -J^ phase 


9^3/16 ^ 0.97 


0.91 



Table 1: Comparison of scalings of boundaries for black rings and plasmarings. 

plasmaring whenever it exists. 

This situation appears to lend itself to a description in terms of a Landau diagram, 
with the entropy given by a function of the (unidentified) order parameter that has two 
maxima (the plasmaball and the large plasmaring) separated by a single minimum (the 
small plasmaring). This analogy suggests - and we conjecture that - the small plasmaring is 
always dynamically unstable, while the plasmaball and large plasmarings are dynamically 
stable with respect to axisymmetric fluctuations. 

An honest verification of our conjecture would require a study of the spectrum of linear 
fluctuations about our plasmaball and plasmaring solutions, an analysis that we have not 
carried out. In this subsection, however, we present some evidence for our conjecture, using 
the 'turning point' stability analysis of [I2] (see [I3] for discussion and references). 

Consider a (not necessarily stable) equilibrium configuration that changes from be- 
ing stable to unstable under continuous variation. The configurations we apply these 
considerations to are plasmarings; according to our conjecture these rings are stable to 
axisymmetric fluctuations when large but become unstable to the same modes when small. 
At the boundary of stability, the matrix of second derivatives of the entropy with respect 
to off shell variations (or 'order parameters') of the configuration under question devel- 
ops a zero eigenvalue. In the neighbourhood of this special point, a small change in the 
thermodynamic potentials of the solution give rise to a large change in the order parame- 
ter along the zero eigenvalue direction (as such a change is entropically inexpensive). As 
argued in [HHTT], this results in a divergent contribution to the second derivative of the 
equilibrium entropy as a function of equilibrium thermodynamic quantities (for instance 
the angular momentum at fixed energy) proportional to the negative inverse of the small 
eigenvalue. 

It follows that a configuration that changes stability has divergent second derivatives 
of entropy with respect to - say - angular momentum. Moreover the sign of this second 
derivative is positive in the 'more stable' phase and negative in the 'less stable' phase. 
Note that the turning point method gives information about the change in the number of 
unstable directions about a solution, but does not yield information about the absolute 
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number of instabilitiesl"] 



The turning point method is useful because it yields information about stability prop- 
erties, with respect to off shell fluctuations, of phases, using information only about on 
shell variations. It is especially useful in the study of nonextensive systems like black 
holes, for which negative specific heats do not necessarily imply dynamical instability 
(note that we're working with the microcanonical ensemble, unlike the grand-canonical 
considerations of [IH])- This method has been used to study the stability of black rings 
in 5 dimensions [ISlIin]; it suggests that small black rings are always unstable, while large 
black rings are more stable in that context. This result corroborates the explicit linear 
fluctuation analysis about the flat space black rings [TO] . 

We now proceed to apply the turning point method to our plasmarings. Define 



In figHS] we have plotted ip against angular momentum at fixed energy for our ring 
solutions. This graph has a single turning point, precisely at the point at which the large 
ring turns into a small ring. The slope of the curve turns from positive (for the large 
ring) to negative (for the small ring) in upon passing through the turning point, consistent 
with our conjecture about the stability properties of plasmarings. More generally, figlTSl 
is qualitatively similar to the equivalent graph of [131 fig-6(b)] for black hole and black 
rings in flat 5 dimensional space, except that the large black ring curves back down as we 
increase L. This difference has no impact on stability analysis, as the turning point method 
links instabilities to vertical tangents rather than horizontal tangents (even though a heat 
capacity/susceptibility changes sign as one crosses a horizontal tangent). 

In conclusion, the turning point method indicates that the small ring has an additional 
instability as compared to the large ring. Note that it is perfectly possible that both the 
large and the small ring are unstable, for example to fluctuations that break rotational 
symmetry. 

6 Four dimensional plasmarings 

In the rest of this paper we turn to a consideration of localised plasma configurations 
in certain massive 4 and 5 dimensional field theories obtained by compactifying related 
5 and 6 dimensional CFTs on a Scherk-Schwarz circle. Although the field theories in 

^^Moreover, this method only Unks vertical tangents - and not vertical asymptotes - in the graph of the 
first derivative of entropy with respect to (say) angular momentum vs. angular momentum (a conjugacy 
diagram) to instabilities, as vertical asymptotes occur at boundaries of equilibrium solution space instead 
of separating solutions of differing degrees of stability. 




(56) 



(57) 
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Figure 13: Conjugacy diagram: ip as a function of angular momentum, L for fixed energy, 
E = 40. 



question are not gauge theories (e.g. the 5 dimensional massive theory could be obtained 
by compactifying the (0,2) theory on the world volume of an M5 brane on a Scherk-Schwarz 
circle) they undergo first order 'deconfining' transitions and the high temperature phase 
of these theories admits a fluid dynamical description. The fluid conflgurations we will 
construct are dual to localised black holes and black rings in Scherk-Schwarz compactifled 
AdSe and AdSy respectively. 

In this section we study stationary solutions of fluid dynamics in 3+1 dimensional fleld 
theories. Our study will be less thorough than our 3 dimensional analysis above; we flnd 
solutions analogous to those in 3 dimensions, but we postpone the complete parameteri- 
sation and study of the thermodynamic properties of these solutions to future work. In 
Appendix [X] we have derived the equations relevant to stationary fluid flow in 5 dimensions, 
but we leave the study of their solutions (and their higher dimensional counterparts) to 
future work. 



6.1 Stress tensor and equations of motion 

In this section we set up the equations of motion of our fluid. We proceed in direct imitation 
of our analysis of = 3 above. We use the metric 



= -dt2 + dr^ + r^d^^ + dz^ 



(5J 



This gives the same non-zero Christoffel symbols as before ffT^ . We choose the origin so 
that r = is the axis of rotation and there is a reflection symmetry in the plane z = 0. 



For our conflgurations, = 7(1, 0, u, 0) with 7 



:i 



-1/2 



We assume that the 



surface can be described by /(r, z) = z — h{r). In the interior of the fluid, This leads to the 
stress tensor 

/72(p + cuVP) 



" perfect 



V 



iMp + p) 





p 







2i(^Vp + P) 









(59) 
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and the surface stress tensor 



rpfJ.V 

surface 



a6{z — h{r)) 



n- 


h /i'(r)2 








\ 







-1 





-h'{r) 










l+/i'(r)2 





\ 





-h'{r) 








(60) 



Just as in (i = 3, the only potentially nonzero term in ^dissipative proportional to 

^ . As in = 3, it will turn out that this quantity vanishes on our solutions, so we 

simply proceed setting T^-,,,p^uve to zero. 

The equations of motion, V^T'^'^ = 0, reduce to 



_d_ 

dr 







dP 



2 

uj r 



dr 1 — cj^r^ 
dP 



{p + P)T 2aHh'{r) d{z - h{r)) 



(61) 



= — ±2aH5{z-h{r)), 

where the upper sign refers to the upper {z > 0) surface and 

rh" + h'(l + h'^) 



H = T- 



2r(l + /i'2)3/2 



(62) 



is the mean curvature of the surface [20] . 
6.2 Solutions 

Our analysis so far has been rather general; to proceed further we use the equations of 
state ([111). We define dimensionless variables as before 



auj 



UJ = — , V = ujr ^ h{v) = ujh{r) . 
Po 



(63) 



Using the equation of state ( jT4j) . we can rewrite (!6T|) in the bulkinterior of the fluid as 

1 dp 5v 

~ ' ' ' (64) 



p — Podv 1 — t>2 

^ (p(.)-Po) {l-v'f 

where K is an integration constant. This means that the pressure and temperature are 



P = Po 



K 



(1 _ ^;2)5/2 

(notice that this justifies our neglect of heat flow). 



T = 7 



a 



(65) 
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Integrating (!UT|) across an outer surface gives 

P = 2aH or — r - 1 = -u ^ — . 66 

{l-v^y/^ + /l'2)3/2 ^ ' 

This can be integrated once to give 

vh' K C 

+ — + ^ , (67) 



where C is another integration constant. 

If we introduce a parameter / equal to the distance along the surface, measured from 
(w, h) = {vo, 0), we have ^ = — \/l + h'"^. Then ( 1U7|) can be written as 

dh 2K-3(t;2 + 2C)(l -t;2)3/2 



d/ Q^v{l - 1;2)3/2 



dt; y^36S32^;2(l _ ^;2)3 _ [2K - 3(t;2 + 2C)(1 - i;2)3/2]^ 



(6J 



d/ 625t;(l - 1;2)3/2 

dh _ 2K - 3{v^ + 2C)(1 - t;2)3/2 

y^36C32i;2(i - y2)3 _ [2i^ _ 3(^2 + 2C)(1 - ?;2)3/2]2 

It follows that the outer surface of our plasma configuration is given by 

/ , . 2i^-3(x2 + 2C)(l-x2)3/2 



\ y'36S32x2(l-x2)3-[2i^-3(a:2 + 2C)(l-x2)3/2]^ 

Of course this only makes sense provided 

6i5x(l-x2)3/2 > |2A^-3(x2 + 2C)(l-x2)3/2| \/ x e {v,Vo). (70) 

Note also the conditions p > Po =^ K > and, of course, < f i < Vo- 

Inner boundaries to the plasma configuration (if they exist) obey the equation P = 
—2aH. The equivalent of (l68ll . with a new integration constant D replacing C (the inte- 
gration constant i^' is a property of the plasma, not the surfaces), is 

dh _ _2K - 3(^2 + 2D){1 - i;2)3/2 
dT ~ 625w(l - 1;2)3/2 ' 



dt; y'36252t;2(l _ ^,2)3 _ \2K - 3(t;2 + 2D){1 - i;2)3/2]2 

dT " 6i5t;(l - i;2)3/2 

d^ 2i^^-3(i;' + 2D)(l- 17^)3/2 

^36S52t;2(i _ ^;2)3 _ [2K - 3(i;2 + 2D){1 - v'^f/^f 



(71) 
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Horizon topology 


Plasma topology 


Object 


^4 

X 52 
52 X 51 X 


53 

52 X 

51 X 52 

5^ X 51 X 51 

None 


Ball 
Ring 

Hollow ball 
Hollow ring 
None 



Table 2: Topologies of gravity and plasma solutions 

The profiles of such boundaries may be obtained by integrating the equation above. 

Even before doing any analysis, we will find it useful to give names to several easily 
visualised, topologically distinct fluid configurations. 

Ordinary ball: v'{l) = h{l) = Q Sit v = Vo- h'{l) > for < < Vo- h'{l) = at = 0. 

Pinched ball: v'{l) = h{l) = at = v^. h'{l) > for < < ^m- ^'(0 = at = 
h'{l) < for < < t;^. h'{l) = at = 00 

Ring: v'{l) = h{l) = at i; = v^. h'{l) > for i;^ <v< v^. h'{l) = at w = y^. h'{l) < 
for vi < V < ih^. v'{l) = h{l) = at = fi, where Vi < ih^ < Vo- 

Examples of these surfaces can be seen in figs JTMTTl Each of these solutions could have 
lumps of fluid eaten out of them. We will use the terms 

Hollow ball: A ball (pinched or ordinary) with a ball cut out from its inside. 

Hollow ring: A ring with a ring cut out from its inside. 

Toroidally hollowed ball: A ball with a ring cut out from its inside. 

It is easy to work out the horizon topology of the gravitational solutions dual to the 
plasma topologies listed above. [22] have obtained a restriction on the topologies of horizons 
of stationary black holes in any theory of gravity that obeys the dominant energy condition; 
any product of spheres obeys the conditions from their analysis. Although the dominant 
energy condition is violated in AdS space, in table [21 we have listed all 4 dimensional 
horizons that are topologically products of lower dimensional spheres, and note that all 
but one of these configurations is obtained from the dual to plasma objects named above 
{B^ is a ball, B^ is a disc and 5^ is an interval). The last one, T'^ = x x x , is 
a marginal case of the theorem. 

In the rest of this section we will determine all stationary, rigidly spinning solutions of 
the equations of fluid dynamics described above. 

^^Black holes with wavy horizons in six dimensions and above were predicted in [21j . 
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6.2.1 Ordinary ball 



We search for solutions of (|69l) for which h'{v) vanishes at w = and blows up at the 
outermost point of the surface vo', we also require that h decrease monotonically from to 
Vo- The first condition sets K = 3C. The condition that v'{l) is zero at Vo may be used to 
determine u; as a function of Vo and K from the linear equation 

2K - {3vl + 2K){1 - vlf'^ = 6uv,{l - vlf' , (72) 

(the choice of positive square root comes from the fact that the LHS above is positive). 
Note that the numerator of the formula for h'(v) be written as 



2h n .^m(K ^'^'i^-^'f" \ 



and 2 (l — (1 — f 2)^/2^ > 3^12(1 — t)2)3/2^ Thus, K > 1 guarantees our monotonicity re- 
quirement. From (1651) . we see that this also ensures that the pressure is positive throughout 
the ball. 

In summary, the full set of ordinary ball solution is obtained by substituting C = K/3 
and u = Lj{K,Vo) (obtained by solving fl72|) ) into fl69|) . This procedure gives us a ball 
solution for any choice oi K > 1 and > 0. 

In figs JTUllTTI we present a plot of the profile h{v) for the ball solution at Vq = 0.8, 
K = 1.5. 



6.2.2 Pinched ball 

The pinched ball satisfies all the conditions of the ordinary ball except for the monotonicity 
requirement on h{v); in fact the function h{v) is required to first increase and then decrease 
as V runs from to fo. It follows that C and uj for these solutions are determined as in the 
previous subsubsection (C = K/3 and u from ( I72l) ) however the requirement h"{v) > at 
V = forces K < 1. This ensures that h'{v) > at small v and h'{v) < at larger v. It 
also ensures that the solution has negative pressure at the origin and positive pressure at 
the outermost radius. 

Not every choice of {K,Vo) G [0, 1], however, yields an acceptable pinched ball solution. 
As we decrease Vo from 1, at fixed K, it turns out that h{0) decreases, and in fact vanishes 
at a critical value of Vq. Solutions at smaller Vo are unphysical. The physical domain,in 
{K, Vo) space is given by the inequality 

2K - {3v^ + 2K){1- i;2)3/2 ^ ^ 

36u^v^{l - v^ f - [2K - (3t'2 + 2K){1 - t;2)3/2]2 

(73) 

We should also ensure that (1701) is not violated, i.e. 

Q{vo,K)= ini \3Qu\>\l-vy-[2K-{3v^ + 2K){l-vy/^f}>0. (74) 
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The boundary of the domain permitted by (I7S1) is plotted in figHH We have also plotted the 
boundary of the region where fl74l) is violated. We see that (1731) is the stricter constraint, 
and that (|7^ is not violated for ordinary balls either. The full set of pinched ball solutions 
is parameterised by values of Vo and K in the region indicated in figHH 
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Figure 14: Domain of ball solutions. 



In figs IT6|T7l we present an example of the profile h{v) for the pinched ball solution at 
parameters Vo = 0.8, K = 0.55. 

6.2.3 Ring 

The plasma of the ring configuration excludes the region v < vi; as this region omits v = 0, 
K and C are not constrained as before. 

As v'{l) vanishes at Wi,fo we have the following constraints 

2K - 3{v^ + 2C)(1 - v^f' = -Quv,{l - v^ f" , 

2K-?,{vl + 2C){l-vlfl' = Q::^Vo{l-vlfl\ ^ ' 

the choice of negative /positive square roots comes from the requirements that h'{l) < 
a.t V = Vi and h'{l) > at v = Vo- These equations may be used to solve for C and uj 
as a function of K, vi, Vo- K{vo, vi) may then be determined from the requirement that 
h{vi) = h{vo) = 0, i.e. 

dh ^ r 2J^-3(t;2 + 2C)(l -t;2)3/2 

— dv = - I \ )\ ) dv = Q. (76) 



dv 



y^3625V(l - ^2)3 _ [2A' - 3(t;2 + 2C)(1 - i;2)3/2] 



In practice, it is easier to first eliminate K and C using (1751) . then substitute v\ = cDf;, 
Vo = ojTo and use (1761) to solve for uj at fixed f\ and fo- after this, one can determine v\ 
and Vo from tD, fi and Tq. We present a 3 dimensional plot of K as a function of v^ and Vo 
for 1 < fo < 10, 0.1 < fi/fo < 0.9 in figHSl All of these solutions have > 0, as required. 
Unfortunately, with this method, one cannot see if there is a physically acceptable solution 
for the whole range of < t"! < fo < 1. It appears that there is a solution for every value 
of fi < fo- 

In figs JT6|T71 we plot the profile function h[v) for the ring solution at parameters fi = 10, 
fo = 20. 
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Figure 15: K as a function of V[ and for ring solutions. 



6.2.4 Hollow ball 

In this subsection we will demonstrate the non-existence of rigidly rotating hollow ball 
solutions to the equations of fluid dynamics. Let us suppose such a solution did exist. The 
inner surface must have vanishing gradient at f = 0; this sets D = K/3. Now let the 
outermost point of the eaten out region he v = Vq. The inner surface must have a vertical 
tangent at Vo- This also implies that the outer surface also has a vertical tangent at Vq 
(the condition for a vertical tangent is identical for an outer or inner surface). However, 
such points saturate the inequality (1701) and, as discussed in §§ §6.2.21 this never happens 
in the interior of a ball. It follows that hollow ball solutions do not exist. 

6.2.5 Hollow ring and toroidally hollowed ball 

Let us first consider the possibility of the existence of a toridally hollowed ball solution. 
Let the innermost and outermost part of the hollowed out region occur aXv = v\ and v = Vo 
respectively. Let us define a{v) = 6ujv{1 — u^)^/^ and b{v) = —2K + 3(f^ + 2D){1 — f ^)^/^ 
where D is the integration constant for the hollow. From fl71l) it must be that 

a{vo) = b{vo) a{vi) = -b{vi) \b{v)\ < \a{v)\ Wv e {ih, Vo) 

For these conditions to apply, b{v) must start out negative at v = Vi, increase, turn positive, 
and cut the a{v) curve from below at f = Vo- We have performed a rough numerical scan 
of allowed values of parameters {K,u},D); it appears that this behaviour never occurs 
(although we do not, however, have a rigourous proof for this claim). For all physically 
acceptable values of parameters, the curve b{v) appears to either stay entirely below a{v) 
or to cut it from above 1^ 

^•^We emphasise that this behaviour appears to be true only for oj > ujminiK) where ojmin{K) is the 
smallest allowed value of ui at fixed K (see figlM]). It is easy to arrange for b{v) to cut a{v) from below 
when uj is taken to be arbitrarily small at fixed K and D, but this is unphysical. 



32 



(J 4 U.U U.O VJ.J-'i U.J-U U.J-O U.iVI M.ii U-i-i 



Figure 16: Profile of tlie surface of an ordinary ball, pinched ball and ring. 




Figure 17: 3D plot of the surface of an ordinary ball, pinched ball and ring. 



These considerations, which could presumably be converted into a proof with enough 
effort, lead us to believe that the existence of hollow balls is highly unlikely. We believe that 
similar reasoning is likely to rule out the existence of hollow rings, although this is more 
difficult to explicitly verify, as our understanding of the parameter ranges for acceptable 
ring solutions is incomplete. 

In order to understand intuitively why hollow rings and toroidally hollow balls are 
unlikely, note that the pressure at the inner and outermost parts of the hollowed out 
region is given by 

P{v,) = pou (^-KUJ + , P{vo) = Po^ - ^) ' 

where v"^fj. is positive and v"^^;,^ is negative. Provided that | and are not 

drastically different, we would require P{vi) > P{vo)- However, the pressure increases 
monotonically with radius. 

In conclusion, we strongly suspect, but have not yet fully proved, that the full set 
of rigidly rotating solutions to the equations of fluid dynamics in ci = 4 is exhausted by 
ordinary balls, pinched balls and rings. 
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7 Discussion 



In this paper we have emphasised that the AdS/CFT correspondence imphes a duahty 
between nonsingular classical gravitational solutions with horizons, and solutions to the 
boundary equations of fluid dynamics. This connection has previously been utilised by 
several authors to obtain gravitational predictions for various fluid viscosities and conduc- 
tivities (see, for instance, [3] and references therein). The new element in our work is the 
incorporation of boundaries separating the fluid from the vacuum into the Navier-Stokes 
equations. This feature (which relies on the explicit gravitational construction of the do- 
main wall in [2]) allowed us to study stationary finite energy lumps of plasma, which are 
dual to localised black holes and black rings in the bulk. 

All our work (and easily imagined generalisations) apply to confining field theories. 
Stationary black holes in such backgrounds sit at the IR ends of the geometry; the boundary 
shadow of such black holes is a lump of deconfined fluid of size R + 0{A~J^p). The fluid 
dynamic equations accurately describe such lumps only when R ^ ^gap; the same 
limit the fluid yields an approximately local representation of the horizon. Consequently, 
the AdS/CFT correspondence provides an approximately local fluid description of horizon 
dynamics in the long wavelength limit. This result is strongly reminiscent of the Membrane 
paradigm of black hole physics [5l[23l[2l], and may constitute the precise version of this 
idea in the context of asymptotically AdS spaces. 

All the specific results of this paper are based on the equations of state ffT^ . which 
are valid for the high temperature phase of Scherk-Schwarz compactified conformal field 
theories (dual to gravity in Scherk-Schwarz compactified AdS space). However the only 
qualitative feature of this equation of state that was important for the existence of the 
solutions of this paper is that the fluid pressure vanishes at finite energy density. In 
fig|9] we have plotted the maximum and minimum values of lnT/7^ in our solutions, as 
a function of energy and angular momentum. Note that at large values of charges (and 
away from extremality bounds) each of these quantities tends to zero. This demonstrates 
that over large classes of our solutions, the fluid temperature is always close to the phase 
transition temperature. As a consequence such solutions 'sample' only the fluid equation 
of state only in the neighbourhood around the zero pressure point, and so would exist in 
any fluid whose pressure vanishes at finite energy density. 

Our results suggest several directions for future research. It would be interesting to 
analyse the stability of small fluctuations about the solutions presented in this paper. As 
we have mentioned in § §5. 51 we expect the small ring to be unstable to axisymmetric fluid 
fluctuations, while we expect the ball and the large ring to be stable to such fluctuations. 
However, it is quite possible that such an analysis would reveal that the large ring solutions 
of this paper have a Plateau-Rayleigh type instability that maps to Gregory-Laflamme 
instabilities (see also [1])0 of the dual bulk solutions. 

Although we have not mentioned this in the text, there exists a scaling limit in which 

^"'We thank T. Wiseman for suggesting this. 
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the thin plasmarings solutions simphfy greatlyo In this hmit (C; — with fixed), the 
3D plasmaring reduces to a straight strip of moving fiuid. The fiuid pressure vanishes 
on this strip, and the fiuid velocity is constant across the strip (more precisely Vo = 
Vi + ^^;^u} + 0{uj'^) SO that To — fj = ^^^ + 0{lu)). Similarly, there should exist scaling limit 
under which the 4 dimensional plasmaring reduces to an infinite stationary cylinder, with 
fiuid fiow along the axis. Various dynamical properties of large rings (e.g. the potential 
Gregory-Lafiamme type instability alluded to in the previous paragraph) will probably 
prove easiest to study in this scaling limit. 

It should also be relatively straightforward, and rather interesting, to more fully anal- 
yse the thermodynamics of the four dimensional solutions presented in this paper. This 
thermodynamics may have interesting features; for example, it has been suggested that 
ultra-spinning black holes in six dimensional fiat space are unstable [21] and it would be 
interesting to see if the same is true of our (pinched?) plasmaballs. 

An extension of our work to obtain the moduli space of five and higher dimensional fiuid 
configurations - and so seven and higher dimensional gravitational black solutions should 
also be possible (though analytic solutions may be harder to obtain in higher dimensions). 
Such an extension would yield interesting information about horizon topologies in higher 
dimensional gravitational theories. An obvious conjecture based on intuition from fiuid 
fiows would be that the full set of stationary fiuid solutions in five dimensions appear in 
three distinct topological classes; solutions whose bulk dual topologies would be 5"^, x 
and X X S^. The reason one might expect the last solution is that in five (but no 
lower) dimensions, it is possible to have solutions that rotate about two independent axes; 
these two rotations should be able to create their own distinct centrifugal 'holes', resulting 
in the above topology. It would be very interesting to check whether this conjecture is 
borne out. 
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Appendices 



A Five dimensional plasmarings 

In the bulk of this paper we have presented an analysis of stationary fluid configurations 
of the three and four dimensional fluid flows. The analysis of analogous configurations in 
one higher dimension has an interesting new element. The rotation group in four spatial 
dimensions, 5*0(4), has rank 2. Consequently a rotating lump of fluid in five dimensions 
will be characterised by three rather than two conserved charges (two angular momenta 
plus energy). When one of the two angular momenta is set to zero, it seems likely that 
the set of stationary solutions will be similar to those of the four dimensional fluid; in this 
limit we expect ball and ring configurations whose dual bulk horizon topologies are 
and X respectively. However turning on the second angular momentum on the ring 
solution could centrifugally repel the fluid away from the second rotational axis, leading 
to a fluid configuration with dual bulk horizon topology x x . Such configurations 
have not yet been discovered in gravity, and it would be exciting to either construct them 
in fluid mechanics, or to rule out their existence. 

In this appendix we set up and partially solve the equations of stationary fluid flow in 
five dimensions. While the stationary equations of fluid dynamics are trivial to solve in the 
bulk in every dimension, boundary conditions are harder to impose in higher dimensions. 
In the particular case of 5 dimensions, the imposition of these boundary conditions requires 
the solution of a 2nd order ordinary differential equation, that we have not (yet?) been able 
to solve. It may be that a full study of this case would require careful numerical analysis, 
which we leave to future work. In the rest of this appendix we simply set up the relevant 
equations, and comment on the dual bulk interpretations of various possible solutions. 

Consider a fluid propagating in flat five dimensional space 

ds^ = -dt^ + drl + rld(f)l + dr^ + r^d02 • 

Consider a fluid flow with velocities given by = 7(1, 0, cui, 0, 1^2)5 where 7 = (1 — — 
f|)~^/^, vi = uiri and V2 = 102^2- Let the fluid surface be given by /(ri, r2) = r2 — h{ri) = 0. 
The stress tensor evaluated on such a fluid configuration is given by 

7V(p + P) \ 


7V^2(p + ^) 
P 

'^y2P+ii-vi)p)j 



perfect 







\^ '-l'^UJ2{p + P) 



vi)P) 



YMp + p) 

p 

'^{v!p+{l-vi)P) 
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The nontrivial equations of motion that follow from 
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justifying this when we find that T oc 7. 
take the form 








dr^P 



aGh'{n)5{r2 - h{n)) 



(77) 



where 



G 



nhh" + {1 + h'^){hh' - n) 



rih{l + /i'2)3/2 

Using the equation of state (fill) , the equations of motion in the fluid interior 

dp , 2vi (^P n., \ 2f2 



dvi 



Hp- po)- 



are easily solved and we find 



(p - po)(l -v'i- vlf = 5Kpo P = po 



K 



T 



1 



Kpo 



a 



nl/6 



Integrating the equations of motion (|771) across the surface we obtain the condition (the 
upper sign should be used for upper surfaces) 



K 



[1 - ufrf 



a nhh" + {i + h'^){hh' -n 



Po 



ri/i(l + /l'2)3/2 



(75 



Unfortunately we have not yet been able to solve this equation; we postpone further 
analysis of (ITHl) to future work. In the rest of this appendix we qualitatively describe 
possible types of solutions to these equations, and their bulk dual horizon topologies. 

In figim we have sketched some possible topologies for these solutions. The first touches 
both the ri = and r2 = axes and we refer to this as a ball. The second type only touches 
one of these axes and we refer to this as a ring. The third type touches neither of the axes, 
as the plasma has the topology of a solid three-torus we refer to this as a torus. 

Each of these could be pinched near either axis, and there could be hollow versions 
(though the considerations of ^ make it seem unlikely that hollow configurations will 
actually be solutions). 
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Figure 18: Topologies of five dimensional solutions. 



The horizon topology of the dual black object can be found by fibering three circles over 
the shapes in figlTSl One of these circles degenerates at each axis (the angular coordinates 
01 and 02); and the other degenerates on the fluid surface (the Scherk-Schwarz circle). The 
topologies generated are: 
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Object 








Ball 






53 X 51 


Ring 






52 X 51 X 


Torus 


X 




fii X 53 


Hollow ball 


S^xS^x 




X X 


Hollow ring 


S^x S^x S^x 




xS^ xS^ X 


Hollow Torus 


S^x S^x S^x 
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None 
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B Notation 



We work in the ( — h +) signature, /i, ly denote space-time indices. 



Symbol 


Definition 


Symbol 


Definition 




Plasma free energy 


/ 


Free energy density 


S 


Plasma energy 


P 


Proper density 


S 


Plasma entropy 


s 


Proper entropy density 


r 


Plasma temperature 


P 


Pressure 


Po 


Plasma vacuum energy 


% 


Deconfinement temperature 


a 


see dn]) 


Pc 


Deconfinement density 




Stress tensor 




^ = 7(1,^^) 


a 


Surface tension 


7 




fix) 


Surface at /(x) = 




Angular velocity 


f/j. 


d,f/Vdf-df 


V 


ojr 




see (HD 




Bulk, shear viscosity 




Heat flux 


K 


Thermal conductivity 
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Symbol 


Definition 


Symbol 


Definition 


U) 


cruj/po 


f 




To 


Outer radius 






r; 


Inner radius 








Vo/uJ 


9±{v) 


Boundary conditions (|36l) 


fi 


Vi/uJ 






E 


Ball/ring energy 


S 


Ball/ring entropy 


L 


Angular momentum 


T 


Ball/ring temperature 


Vt 


Ball/ring angular velocity 


An 


Vahdity criterion (1541) 


E, L, S*, T, 


See gID 




{dS/dE)i 


0,A,B,C 


Existence regions, see figE] 




(dS/dL)^ 


6,A,B,C 


Existence regions, see figH 




Boundarv scalings in ^^^5.41 


h 


Surface height z = h{r) 


h 


ujh 


H 


Mean curvature 


K, C,D 


Integration constants in [|6] 


Ah,Q 


Consistencv conditions (173(741) 


G 


Curvature in ^JSl 
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